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A global two-neutrino oscillation fit combining Super-Kamiokande solar neutrino 
data with the solar neutrino rates measured by Homestake, Gallex/GNO, SAGE 
and SNO prefers a single allowed area, the Large Mixing Angle solution, at about 
95% confidence level. The mass 2 difference between the two mass eigenstates is 
Am 2 «3-25xfO~ 5 eV 2 , the mixing angle 6 is tan 2 9 R*0.21-0.67. 



1 Introduction 

Neutrino flavor oscillations are currently the most favored way to understand 
the larger solar 8 B neutrino interaction rate measured by Super-Kamiokandad 
(SK) with respect to the 8 B rate reported by the SNO-collaboratiorio. The 
SNO rate includes only e-type 8 B neutrinos while SK measures solar neutri- 
nos via neutrino-electron elastic scattering which has a small sensitivity to 
other (active) flavors as well. The SK excess is then interpreted as a hint for 
appearance of other active flavors in the 8 B neutrino flux (born in the sun 
as purely e-type) . Neutrino oscillations also explain the xlcficit of both rates 
as well as the deficits, of the "Chlorine rate" (HomestakeEl) and the "Gallium 
rates" (Gallex/GNO^, SAGEQ) with respect to the Standard Solar ModeB 
(SSM). Chlorine and Gallium rates include only e-type solar neutrinos. 

The large neutrino flavor mixing between the second and third generation 
inferred from atmospheric neutrino dataO and the absence of an oscillation 
signal in the CHOOZ reactor neutrino experimentB requires the neutrino flavor 
mixing between the first and third generation to be small. Solar neutrino 
oscillations can therefore be approximated by a two-neutrino description with 
the parameters 9 (mixing angle) and Am 2 (difference in mass 2 ). The mixing 
angle ranges between and tt/2, since Am 2 is defined to be positive. For 
Am 2 between « 10 _8 eV 2 and ~ 10 _3 eV 2 , the matter density in the sun and 
earth can strongly affect the oscillation probability. For 8 < 7r/4, resonant 
enhancement of the oscillations (MSW effectEl) can occur in the sun. On 
the "dark side" of the parameter space (6 > ir/4), an anti-resonance can 
suppress the oscillations. Below « 10~ 9 eV 2 (quasi- vacuum/ vacuum region), 
the oscillation probability is more affected by the oscillation phase than by 
matter effects. 



2 Parameter Estimation for Two-Neutrino Oscillations 



The simplest way to constrain the oscillation parameters is a \ 2 fit to all e- 
type neutrino rates, that is Gallex/GNO and SAGE (combined into a single 
"Gallium" number), Homestake and SNO. The dotted lines of figure 0a) show 
the allowed regions-a.t 95% C.L. (which can be understood as an "overlay" of 
the allowed regionst2l from the three experimental rates mentioned above and 
their SSM-based predictions). The area above Am 2 « 10~ 5 eV 2 near maximal 
mixing is called the Large Mixing Angle (LMA) solution. The Small Mixing 
Angle Solution (SMA) is located between Am 2 « l(T 6 eV 2 and « lCT 5 eV 2 
at tan 2 9 ss 10 -3 . The LOW solution is the large Am 2 part of the extensive 
region(s) between ss 10 _7 eV 2 and ~ 10 _9 eV 2 . The lower part of this region is 
called the quasi-vacuum (quasi- VAC) solution. The vacuum solutions (VAC) 
are below w 10~ 10 eV 2 . All regions have similar \ 2 (LOW fits slightly worse). 
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Figure 1. SK Zenith Spectrum/SNO (a) and global fit (b) contours at 95% C.L. The 
dark-grey area in a) is excluded by the shape of the SK Zenith Spectrum. Adding the SNO 
rate measurement results in the (light-grey) allowed areas. Both regions do not depend 
on neutrino flux predictions by the SSM. Overlaid are the allowed regions (inside dotted 
lines) based on the Gallex/GNO, SAGE, Homestake and SNO rates and the SSM. In b) the 
allowed areas from a global fit (adding the rates of Gallex/GNO, SAGE and Homestake to 
the SK Zenith Spectrum-SNO fit) are predominantly LMA solutions. 
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To break the degeneracy in \ 2 of these regions, more information is needed 
than provided by the suppression of the e— type rates. Indeed, neutrino oscil- 
lations not only cause a 'disappearance' of the e— type fluxes but also produce 
other neutrino flavors (appearance). Furthermore, they can induce distortions 
of the neutrino spectrum and time variations of the solar neutrino fluxes. In 
the MSW region, the time variations arise from matter effects inside the earth 
(daily variations). In the vacuum regions, the time variations are the conse- 
quence of the change of the oscillation phase due to the yearly variation of 
the oscillation baseline (distance between sun and earth). 

SK can search for spectral distortions through analysis of the spectrum of 
the recoiling electron as well as daily or year^titne variations. The analysis of 
the shape of the SK "zenith angle spectrum"li3'liil combines spectrum and daily 
variation analyses. Neither spectral distortion nor daily variation was found. 
The centerpiece of this work is the preliminary SK zenith angle spectrum 
based on the entire SK data set (1496 days). The numerical results are listed 
in tables [j] and [2| Figure 0a) shows the excluded regions (dark-grey) from the 
SK zenith angle spectrum (for an explanation of the fit method see Appendix, 
eq. ||) . The SMA and i£AC regions are excluded since they require a distorted 
8 B neutrino spectrumEj. The lower part sii the LMA and the upper part of 
the LOW solution predict daily variationstj and are therefore disfavored. The 
excluded areas are independent of the SSM neutrino fluxes. The SK zenith 
angle spectrum shape is best described by quasi- VAC solutions. 

If combined with the SNO ratefl of 0.346lg;g^|xSSM (Appendix, eq. |), 
the SK rate provides a probe for the appearance of other neutrino flavors. The 
SK rate of 0.465^;^ x SSM (see table |l|) exceeds the e-type rate inferred from 
the SNO measurement by more than 3cr. If this is interpreted as appearance, 
the other flavors contribute about 25% to the SK rate and 70% to the 8 B flux 
(the SK cross section for the other flavors is six to seven times smaller than for 
e-type neutrinos). The two light-grey allowed regions of figure |l|a) are-based 
on the combined fit to the SK zenith angle spectrum and the SNO rateO. One 
region contains the upper part of the LMA solution, the other region contains 
the lower part of the LOW solution and the quasi- VAC solutions (best-fit). 
The allowed areas are still independent of the SSM neutrino fluxes. nn 

The Appendix (eq. 0) explains the method to n add the GalliumQiJ 
(74.8t|;£SNU or 0.5841^ x SSM) and Chlorine ratesQ (2.56 ± 0.23SNU or 
0.337 ± 0.030 x SSM) to the fit. The allowed areas of this global fit (shown in 
figure Qd) look quite different when compared with the e— type rate fit: Only 
the upper part of the LMA survives (and two tiny quasi- VAC solutions). The 
disappearance of the LOW solution is due to predicted daily variations. The 
SK rate requires a larger 8 B flux than the Chlorine rate allows, so almost all 
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Table 1. SK Rate and Uncertainty for eight energy bins. The rates, statistical and system- 
atic uncertainties (of the spectrum shape) in the third column are presented in units of SSM 
expectation. These systematic uncertainties are assumed to be uncorrelated in energy. The 
uncertainties in the fourth (uncertainty of the 8 B neutrino spectrum), fifth (uncertainty of 
the energy scale of 0.64%) and sixth column (uncertainty of the energy resolution of 2.5%) 
arc fully correlated in energy (but uncorrelated with each other). The combined uncertainty 
(last row) is based on the sum of all bins of the zenith angle spectrum. The combined rate 
has an additional systematic uncertainty of osk = !^2 g % ( excm ding S B shape, energy 
scale and resolution), which was added to the uncorrelated systematic uncertainties of the 
spectrum shape. 



Bin 


Range [McV] 


Rate+stat+syst [SSM] 


S B Spectrum 


E-Scale 


E-Resol. 


1 

2 
3 


5.0-5.5 
5.5-6.5 
6.5-8.0 


0.4671±0.0404 jTooill 
0.4580+0.0141 locoes 
0.4729±0.0084±0.0065 


+0.04% 
-0.02 /0 
+0.13 % 
-0.09 /0 
+0.41% 
-0.38 /0 


+0.09% 
-0.01 /0 
+0.20% 
-0.16 /0 
+0.63% 
-0.62 /0 


+0.23% 
-0.21 /0 
+0.19% 
-0.17 /0 
+0.17% 
-0.16 /0 


4 
5 
6 


8.0-9.5 
9.5-11.5 
11.5-13.5 


0.4599+0.0093+0.0063 
0.4627+0.0103+0.0063 
0.4621+0.0168+0.0063 


+0.89% 
-0.85 /o 
+1.7% 
-1.6 /o 
+3.1% 
-2.7 /o 


+1.3% 
-1.3 /0 
+2.5% 
-2.4 /0 
+4.4% 
-4.1 /0 


+0.12% 
-0.12 /0 
+0.16% 
-0.18 /o 
+1.1% 
-1.1 /0 


7 


13.5-16.0 


0.5666±0.0390±0.0078 


+5.1 % 
-4.2 /o 


+7.0% 
-6.4 /o 


+3.2% 
-3.2 /o 


8 


16.0-20.0 


0.5554+0.1458+0.0076 


+7.7% 
-5.6 /o 


+10.6% 
- 9.6 /o 


+8.4% 
-7.9 /o 


Comb. 5.0-20.0 


0.4653±0.0047i§;^|| 


+1.15% 
-1.04 /o 


+1.66% 

-1.58 /o 


+0.33% 
-0.34 /o 



Table 2. Subdivision of bins 2—7 according to the solar zenith angle 9 Z . The range of cos 6 Z 
is given for each bin: cos(? z < is 'Day' and cos0 z > is 'Night' ('Mantle' and 'Core'). 
The rates are given in units of O.OQlxSSM. Only statistical uncertainties are quoted. All 
systematic uncertainties (see table hi) are assumed to be fully correlated in zenith angle. 





Day 






Mantle 






Core 


Bin 


-1-0 


0.00-0.16 


0.16-0.33 


0.33-0.50 


0.50-0.67 


0.67-0.84 


0.84-1 


2 


453+20 


442+53 


379+49 


472+45 


522+45 


503+49 


426+52 


3 


474+12 


530+34 


506+30 


438+26 


478+26 


451+28 


439+31 


4 


448+13 


463+36 


470+33 


462+29 


509+29 


461+32 


451+35 


5 


453+15 


449+40 


502+38 


451+32 


473+32 


477+35 


483+40 


6 


477+25 


509+67 


351+55 


391+49 


498+53 


434+56 


521+64 


7 


511+54 


570+150 


831+167 


694+131 


665+127 


441+118 


469+131 



the quasi- VAC regions disappear. The fit does not depend on the 8 B and hep 
neutrino flux predictions by the SSM, which suffer from the largest uncertain- 
ties. However, it depends on the other SSM neutrino fluxes, in particular the 
7 Bc flux (10% uncertainty) and the neutrino fluxes of the CNO cycle (« 20% 
uncertainty). Those fluxes contributes about 15%( 7 Be) and 6%(CNO) to the 
Chlorine rate and 27%( 7 Be) and 7%(CNO) to the Gallium rate in the SSM. 
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Table 3. Parameters for the best-fit points. The probabilities given in the fourth row are 
based on the difference of x 2 (with respect to the minimum). The five rows below show the 
five independentparts of the fit: the Ax 2 from a fit to the shape of the SK zenith spectrum 
(Appendix, eq. 0) and four interaction rates (deviation probabilities are given in units of 
Gaussian standard deviation a). The last four rows show the values of the minimized fit 
parameters. The 8 B and the hep fluxes are free, the 8 B neutrino spectrum shift as well as the 
SK energy scale and resolution shifts are constrained within the systematic uncertainties. 



Solution 


Large Mixing Quasi- Vacuum Low Am 2 Small Mixing 
Angle (LMA) (Quasi- VAC) (LOW) Angle (SMA) 


Am 2 
tan 2 


6.0xl0" s 4.57xl(r 1,J 5.0X10" 8 4.8x10"° 
0.35 2.1 0.83 0.00044 


X* (45 dof; p x2 [%]) 
A X 2 (2 dof;p Ay2 [%]) 


43.4 (53.9) 48.5 (33.4) 51.2 (24.3) 54.2 (16.2) 
0.0 (100.0) 5.1 (7.9) 7.8 (2.0) 10.8 (0.5) 


A XSK (Pa x 2 W) 
Ga Rate [SNU] 

CI Rate [SNU] 
SK Rate [%SSM] 
SNO Rate [%SSM] 


3.4 ( 1.3a-) 3.1 ( 1.3a) 3.9 ( 1.5a) 5.0 ( 1.7a-) 
73.2 (-0.3a) 69.6 (-1.0a-) 68.2 (-1.3a) 75.1 (+0.1a) 
2.97 (+1.8a) 3.18 (+2. 80-) 3.13 (+2.5a) 2.67 (+0.5a) 
46.4 (-0.1a) 44.7 (-1.4a) 44.9 (-1.2a) 44.1 (-1.9a) 
32.8 (-0.7a) 37.1 (+0.8a) 38.5 (+1.3a) 43.8 (+3.1a) 


4>s B [10<7(cm^s)j 
<t>he P [10 3 /(cm 2 s)] 
8 B Spectrum Shape 
SK E-scale/resol. 


5.62 (+O.60-) 3.71 (-1.7a) 4.04 (-1.2a) 2.71 (-2.9a) 
40 21 8 
-0.3a -0.7a -0.1a -l-O.la 
-0.4a/-0.1a -1.0a/ 0.0a -0.1a/-0.2a +0.1a/-0.4a 



3 Results of the global fit 

Table | compares the four smallest local minima of the x 2 describing the fit. 
The best-fit is located in the upper LMA area. The 8 B flux resulting from this 
fit is somewhat higher than expected by the SSM (5.05^^1 X 10 6 /(cm 2 s)) but 
well within the uncertainty. The hep flux is considerably higher than expected 
by the SSM (9.3 x 10 3 /(cm 2 s)); however, the uncertainty of this prediction is 
thought to be very large. The fit agrees with the SK zenith angle spectrum 
moderately well. Figure ^a) shows the SK zenith angle spectrum and this best 
fit. Figure |2|b) gives a magnified view of the LMA region (in a linear scale of 
tan 2 9) . The la (black) and 3a (light grey) contours are also shown. The fit 
accommodates well the Gallium, SK and SNO rates, however the predicted 
Chlorine rate is about 2a too high. This worsens the otherwise very good 
best- fit x 2 - As a consequence, quasi- VAC solutions, which don't fit well the 
interaction rates in general, appear at 95% C.L. 

The quasi- VAC solution describes the SK zenith angle spectrum slightly 
better than the LMA. The fit is helped by a statistical up-fluctuation of the 
last two bins (see figure ||). Together with a down-shift of energy scale and 8 B 
spectrum, the required oscillation minimum can be 'generated' in the data. 



5 



a) 

5 0.7 
co 

CO 

13 0.6 



•All 

• Day 



0.4 L 



• Night 



0.5 



6 8 10 12 14 16 18MeV 



CO 

co o.f 



a 0.6 
0.4 

0.8 
0.6 
0.4 



DAY T NIGHT 



' 5.5-6.5 MeV 



3.5-8.0 MeV 



rkf. 

DAY NIGHT 



1 9.5-11.5 MeV. 



DAY NIGHT 



; 11.5-13.5MeV 



c) -s 

V 10 r 
> 



E 

< 



DAY jjjjjGHT DAY 

0.5 1 0.5 1 0.5 1 

COS9, 





1 68% C.L 
1 95% C.L. 
99.73%C.L. 




0.2 0.4 0.6 0.8 1 1.2 1.4 

d) tan 2 (0) 

C\J I — ' " " ~~ " "■ ' 

S -7 
.g 10 r 

< 



10 



10 



10 



10 



01 23456789 10 
tan 2 (0)x 10 4 




Figure 2. SK Zenith Spectrum (a) and LMA (b), SMA (c) and LOW/Quasi-VAC (d) 
contours resulting from the global fit. The top panel of a) shows the SK energy spectrum 
during the day (open circles), the night (filled circles) as well as day and night combined 
(stars). Between 5.5 and 16 MeV, zenith-angle distributions arc shown (lower 6 panels of a). 
The dotted line in each of the lower panels is the combined rate in that bin. Error bars reflect 
the statistical uncertainty only. Superimposed is an oscillation prediction (tan 2 9 = 0.34 and 
Am 2 =6.0T0 -S eV 2 ) near the best-fit point. In b)-d) the one (black), two (dark-grey) and 
three a (light-grey) contours of the LMA, SMA and LOW/Quasi-VAC solutions arc shown 
in a linear scale in tan 2 9. The best fits are indicated by the white stars. Superimposed 
is the 95%C.L. excluded area from the SK Zenith Spectrum (inside dotted lines). SMA, 
LOW and Quasi- VAC arc disfavored. 
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The rates, however, do not fit well. The fit struggles to accommodate both 
the Chlorine and SK rate. In this region Chlorine and SK disagree by 3a 
about the 8 B flux (which in the un-oscillated SSM contributes about 76% 
of the Chlorine rate). The resulting best-fit 8 B flux falls 1.7a short of the 
SSM prediction. Even though the quasi- VAC solution has an overall 'C.L. 
threshold' of 8%, it is considerably disfavored, when thus checked in detail. 
The surviving (at 98% C.L.) LOW solution fits the rates about as poorly as 
the quasi- VAC solution. The SK zenith angle spectrum fits somewhat worse 
than either LMA or quasi- VAC. The lack of zenith angle variation in the 
SK data reduces the Am 2 (usually around 10 _7 eV 2 ) and worsens the LOW 
best- fit which is already under pressure from the rates. Figure ^|d) shows an 
enlarged view of the LOW and quasi- VAC region. 

A "smaller-than-small mixing angle" solution appears just left of the SMA 
at about the 3a level (see figure |c). The SMA region is defined by the crossing 
of the Gallium and the Chlorine allowed area (for a given 8 B flux) . It therefore 
fits those two rates very well. The SK zenith angle spectrum, however, fits 
the worst of all solutions: The SK spectral data lack the predicted distortion. 
The SNO rate is more than 3a above the measurement, the SK rate is too 
low (2a). The 8 B flux required is very low (3a). 

The 8 B flux comparisons above refer to the SSM value. A recent preci- 
sion measurement of the cross-section of the 7 Be(p,7) 8 B fusion reaction by 
Junghans et alB implies (see M) a flux of 5.93tg;f| x 10 6 /(cm 2 s). If the fit is 
confined to this flux, then all solutions other than LMA are further disfavored. 




Figure 3. SK Spectral distortion for two quasi- VAC solutions (solid line: best fit at 4.57 ■ 
10~ 10 cV 2 , dashed line: 6.68- 10 _10 eV 2 ) and data (open circles). Although the data do not 
support spectral distortions, a shift (stars) in energy scale (— l.Ocr) and 8 B shape {—0.7a) 
and systematic uncertainty of the combined rate ( — 1.4<r) can accommodate the prediction. 
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4 Conclusion 



A global fit to all solar neutrino data was performed using a two-neutrino 
oscillation model. At 95% C.L. only the upper LMA solution and two tiny 
quasi- VAC regions are still allowed. The quasi- VAC regions are disfavored at 
more than 92% C.L. caused by a disagreement about the 8 B flux between SK 
and the Chlorine rate of 3cr; the resulting quasi- VAC 8 B flux fit is l.7a smaller 
than the Standard Solar Model. The LMA solutions above Am 2 = 3- lCT 5 eV 2 
are therefore the only viable solution at 95% C.L. 
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Appendix: Fitting Method 

To estimate and limit oscillation parameters, a x 2 describing the shape of 
the SK zenith spectrum is extended to take into account the measured neu- 
trino interaction-rates by various experiments. Using oscillation probabilities 
(obtained as in 13) and SSM neutrino fluxes, the expected interaction rates 
B°"° (due to 8 B neutrinos) and H° B ° (due to hep neutrinos) in energy bin i and 
zenith-angle bin z are calculated. These rates as well as the SK measurements 
Di z are then normalized by the SSM expectations without oscillations: 

Dobc TTosc i-\ 

7 . i,z 7 i.z 7 _ Ui,z \ 

*i 2 dssm i rxsSM ' v / 



b ^ DSSM i rrSSM ' b ^ DSSM _r_ ZTSSM ' "e,^ DSSM i TT 

i.z ' i,z i,z i,z i.z ' 

The zenith components of the flux difference vector 

A;, z (/3,?7) = (/3 • b ijZ + r] ■ h itZ ) x f(Ei,d B ,Ss,S R ) - di, z 

allow for arbitrary neutrino fluxes (through the free parameters (3 and rf). The 
combined rate prediction is modified by the energy-shape factors 

f(Ei,5 B ,6s,S a ) = fB{Ei,8 B ) x fs(Ei,6 s ) x f R {E u 5 R ) 

with 5b describing the 8 B neutrino spectrum shape uncertainty, 5s describ- 
ing the uncertainty of the SK energy scale (0.64%) and 5 R describing the 
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uncertainty of the SK energy resolution (2.5%) (The shapes fs,fs,fR are 
given in table [l]). All three uncertainties affect the bins of the SK zenith 
angle spectrum in a correlated way. The 7x7 matrices Vi describe statisti- 
cal and energy-uncorrelated uncertainties; the latter are assumed to be fully 
correlated in zenith angle. For any given parameters 5k, the x 2 



xl 
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with c = y:( b s- v f[- b 4 H-Yfl'h) and ? = 



can be written as a quadratic form of 4> an d the curvature matrix Co- The 
minimum is 

8 

Xo,m = X! di ' ^ 1 ' di ~ C °- m With C ° m = & 0,m C Q 0,m" 
i=l 

If the minimum flux vector is scaled by a ( <^> — a x <j) ) then x 2 . constrains 
a to be a = 1 ± \/l/Co,ro- To take into account the systematic uncertainty 
of the SK combined rate osk =_2 6 ^° ( wmcn is completely correlated in 
zenith- angle and energy), Xo is modified to 

Xi=Xo,m+(V-?o,m)Ca (7-~$o, m ) with Ot = ^ ^gff 2 xC . 

V / V / Go,m + *-/&SK 



Xi constrains a to be a = 1 ± y/l/Co, m + that is, the minimum is 

unchanged, but the allowed parameter range for (3 and rj is larger. The total 
X 2 for the SK zenith spectrum shape is then 



XL = Min (x?09, „, fa, f fl) fa) + (^) 2 + (|) 2 + (^) 2 ) (2) 

where all as well as /?, 77 are minimized. n 

A combined fit with the rate measured by the SNO collaboration 2 ] can 
be done by defining b SNO ,h SNO and A SNO (/3,?7) in a similar way as in (Q) and 
form 

$b\ 2 . ( 5 S \ 2 . ( Sr\ 2 , /A SNO (/3,^y) , 



^,»o = Min U ; + (- + - + - + l-^l (3) 



9 



To add the_radio-chemical rate measurements of Homestakei, Gallex/GNoi, 
and SAGEEl (all "Gallium" rates are combined into a single rate), the \ 2 

takes into account the correlations between Gallium and Chlorine mea- 
surements. The 8 B and hep fluxes are constrained by the minimization of 
Xsk-sno) n °t by the Standard Solar Model. The total \ 2 is then a simple 
addition 

X 2 = xLsno + Min ar) min ) + (^^) J ( 4 ) 

where = C SK . SN o, m = ^ SK _ SBOim • Csk-sno • "^sk-sno,™- Numerically, a a is 
found to be about 2.7% (The accuracy of SK and SNO combined) 
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